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ABSTRACT 

Thia  report  contains  the  description  of  a  solution  of  the 
problem  of  self-acting,  gas-lubricated,  partial-arc,  finite-length  Journal 
bearings  for  steady-state  conditions. 

The  results  can  be  used  for  bearings  of  many  types  such  as 
simple,  fixed,  partial-arc  Journal  bearings  (tapered  land;,  axial-grooved 
Journal  bearings,  pivot ed-pad  Journal  bearings,  and  similar  geometric 
shapes. 

Sample  results  are  plotted,  and  their  utilization  for  various 
bearing  types  is  illustrated  by  examples. 
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NOMENCLATURE 

i 

A 

coefficient 

aa 

e 

iteration  increment 

B 

coefficient 

bb 

e 

total  growth  indicator 

c* 

coefficient 

cc 

coefficient 

CL 

w 

load  coefficient  -  p  ^ 

a 

CLT 

total  load  coefficient 

. 

C 

ground  olrara™.  -  Rshoe  - 

C' 

pivot  circle  clearance  -  Rpivot  -  R5haft 

* 

D 

coefficient 

Eij 

Reynolds  difference  equation  coefficient 

e 

shaft  eccentricity  with  respect  to  shoe  circle 

e' 

shaft  eccentricity  with  respect  to  pivot  circle 

F,  . 
ij 

Reynolds  difference  equation  coefficient 

F 

c 

frictional  force 

H 

dimensionless  local  film  thickness  =  k 

HP 

hp 

dimensionless  pivot  film  thickness  “jr 

K 

dimensionless  trailing  edge  film  thickness  =  q 

i 

grid  point  in  9  direction 

j 

grid  point  in  axial  direction 

V 

k 

iteration  number 

L 

bearing  axial  width 
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Nomenclature  (Cont'd) 

1,1  iteration  number 

M  number  of  grid  cells  in  0  direction 

N  number  of  grid  cells  in  axial  direction 

P  dimensionless  pressure  =  P/P& 

P&  ambient  pressure 

R  shaft  radius 

5  coefficient 

T  coefficient 

W  total  load 

Wc  load  component  parallel  to  line  of  centers 

load  component  normal  to  line  of  centers 

Y  axial  coordinate 

or  shoe  angular  extent 

0  angle  between  pivot  and  vertical  axis 

V  relaxation  factor 

6  attitude  angle 

*  eccentricity  ratio  =  e/C 

«'  eccentricity  ratio  =  e'/C' 

Tl  dimensionless  axial  coordinate  =  y/L 

0  angular  coordinate 

2 

A  bearing  number  =  -  ^ 

P  C* 

a 

Ap  bearing  number  based  on  hp,  Ap  = 

HP 
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NOMENCLATURE  (Ooncl.) 
bearing  number  based  on  h^,,  A^,  = 

“t 

absolute  viscosity  of  lubricant 

lead  angle  between  line  of  centers  and  shoe  leading  edge 

same  as  above  for  first  pad  of  grooved  bearing 

angle  between  load  line  and  leading  edge  of  first  pad  of 
grooved  bearing 

angle  between  shoe  leading  edge  and  pivot  point 
shaft  angular  velocity. 
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1.  INTRODUCTION 

We  have  witnessed  how  the  increasing  demands  of  space-age 
technology  have  fostered  impressive  developments  in  the  field  of  com¬ 
pressible  fluid-film  lubrication.  It  is  now  recognized  that  many 
properties  of  gas-lubricated  bearings  make  them  attractive  for  certain 
particular  applications.  Most  of  the  very  simple  bearing  geometries  have 
been  investigated  to  the  extent  that  reliable  design  information  exists, 
both  of  analytical  and  experimental  nature.  However,  it  is  also  true 
that  gas-lubricated  bearings  require  accurate  alignment  of  the  shaft  and 
bearing  assembly  and  Mat  many  of  them  are  prone  to  develop  a  self-excited 
film  instability. 

The  alignment  requirement  has  more  implications  than  simply 
the  need  for  initial  parallelism  of  shaft  and  bearings.  Subsequent 
distortion  of  the  bearing  housing  and  supporting  structures  due  to 
transient  stresses  and  thermal  gradients  may  produce  this  critical 
condition.  Maintaining  a  high  degree  of  alignment  with  some  bearing- 
structure  designs  is  very  difficult  and  may  in  some  extreme  cases  be 
impossible.  .  „  . 

The  question  of  film  instability  can  also  be  very  serious.  It 
has  frequently  been  the  major  obstacle  to  successful  bearing  performance. 

Th6  instability  phenomenon,  especially  of  the  so-called  half -frequency 
variety,  is  most  pronounced  in  the  simplest  geometries  such  as  the  360° 
journal  bearing. 

Half -frequency  whirl  may  be  described  as  a  self-excited  vibration 
caused  by  the  nonlinear  hydrodynamic  bearing  forces.  At  certain  threshold 
speeds,  this  self-excited  instability  will  cause  the  rotor  center  to 
whirl  or  process  in  the  same  direction  as  shaft  rotation  at  a  rate 
approximately  equal  to  £  the  total  rotor  angular  velocity.  This  instability 
is  of  an  altogether  different  nature  from  that  of  the  resonance  encountered 
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at  a  rotor  "critical  speed",  which  is  excited  by  rotor  unbalance.  Rotor 
speeds  above  the  initial  onset  of  whirl  will  result  in  rapidly  increasing 
amplitudes  of  the  whirl  orbits  leading  eventually  to  bearing  failure. 

Thus,  a  properly  designed  gas  bearing  should  have  the  half -frequency 
whirl  threshold  above  the  expected  operating  speed  range. 

Information  now  available  based  on  limited  theoretical  knowledge 
and  scattered  experimental  data,  indicates  that  deviations  from  the  basic 
plain  cylindrical  geometry  improve  the  bearing  whirl  stability  character¬ 
istics  . 

This  report  is  concerned  with  the  determination  of  basic 
characteristics  of  partial-arc,  cylindrical  bearings  which  are  somewhat 
resistant  to  the  difficulties  described  above.  This  in  turn  leads  to  the 
design  of  bearings  with  three  important  geometries:  tilting  pad,  axial 
grooved,  and  fixed  partial-arc  (tapered  land)  journal  bearings.  All  of 
these  geometries  present  favorable  features  with  respect  to  the  stability 
problem,  and  the  pivot ed-pad  bearing,  in  addition,  can  operate  with 
considerable  misalignment.  Theoretical  solutions  of  the  complete  steady- 
state  characteristics  of  gas-lubricated,  partial-arc  bearings  are  not 
presently  available,  and  it  is  felt  that  this  will  make  a  contribution 
to  the  advancement  of  the  state  of  the  art. 

The  film  pressure  generated  by  a  self-acting,  gas-lubricated 
bearing  is  governed  by  the  well-known  isothermal  Reynolds  lubrication 
equation.  For  compressible  lubricants  this  equation  is  a  non-linear, 
non-homogeneous,  partial  differential  equation  with  variable  coefficients, 
for  which  analytical  solutions  are  available  only  in  a  few  limiting 
cases. 

Many  approximate  methods  are  available  for  the  solution  of 
Reynolds'  equation.  One  group  of  these  methods  takes  advantage  of  the 
smallness  of  some  bearing  parameter.  For  example,  a  "small  c  expansion" 
has  been  successfully  attempted  by  the  authors  and  furnishes  data  for  • 
the  region  e  <  0.4.  (Note:  The  symbol  e  is  the  eccentricity  ratio  of 
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thte  journal  in  ita  bearing).  This  work  has  not  been  reported  yet  since 
the  small  «  region  covered  by  the  data  is  not  of  major  interest  at  the 
moment. 

Solutions  have  also  been  attempted  on  the  basis  of  a  "small  A 
expansion",  and  their  validity  region  would  be  restricted  to  cases  where 
only  a  modest  role  is  played  by  lubricant  compressibility.  (Note:  The 
symbol  A  refers  to  the  compressibility  bearing  parameter).  Snell  used 

w*. 

this  approach  but  further  simplified  the  problem  by  taking  the  bearing 

1* 

clearance  to  vary  linearly  from  the  leading  to  the  trailing  edge. 

This  assumption  was  motivated  by  the  fact  that,  for  cylindrical  geometries, 
even  the  "small  A  approximation"  to  Reynolds'  equation  presents  the 
difficult  problem  of  the  solution  of  Hill's  equation. 

Other  methods  of  solution,  such  as  those  utilizing  truncated 
series  expansions  and  integral  fits  such  as  Galerkin's  technique,  would 
be  very  useful  in  obtaining  close  approximations  to  the  load-carrying 
capacity  of  a  bearing  with  any  geometry.  However,  the  basic  drawback  of 
all  these  approximate  methods  is  that,  although  the  overall  pressure  load 
can  be  predicted  quite  accurately,  the  actual  shape  of  the  pressure  pro¬ 
file,  and  consequently  the  position  of  the  point  of  application  of  the 
load,  is  difficult  to  obtain.  This  argument  must  be  considered  quite 
seriously  in  cases  where  the  static  data  have  to  be  used  in  dynamical 
analysis.  Then  derivatives  of  the  pressure  profiles  are  necessary,  and 
the  accuracy  of  the  calculations  depends  considerably  on  the  smoothness 
of  the  data  and  on  the  close  similarity  of  the  approximate  and  exact 
pressure  profiles.  An  additional  difficulty  in  the  application  of  methods 
such  as  Galerkin's  is  encountered  in  the  evaluation  of  the  expansion 
coefficients.  Indeed,  the  non-linearity  of  Reynolds'  equation  makes  it 
necessary  to  employ  numerical  techniques  in  the  integration  of  the  determining 

equations . 

- 

Superscript  numbers  refer  to  the  list  of  references  at  the  end  of  this 
report. 
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On  the  basis  of  these  considerations  the  direct  numerical  method 
of  solution  of  the  exact  equation  governing  the  problem  was  chosen  as  the 
one  best  satisfying  all  of  the  accuracy  requirements.  An  approach  for 
the  use  of  this  analysis  in  the  study  of  -she  dynamical  behavior  of  a 
tilting-pad  bearing  arrangement  is  now  in  progress  at  The  Franklin  Institute. 

In  the  reported  analysis,  the  Reynolds  differential  equation 
is  reduced  to  a  finite  difference  equation  by  standard  methods.  The 
work  follows  the  same  approach  as  the  one  used  by  Gross.  Plots  of  the 
results  are  presented  in  a  form  convenient  for  use  in  design. 
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2.  ANALYSIS 

The  Reynolds  Equation 

The  basic  theory  utilizes  the  well  discussed  assumptions  con¬ 
cerning  neglect  of  the  inertia,  forces  and  temperature  variations  through¬ 
out  the  film.  Then,  the  Reynolds  equation  can  be  represented  in  the 
following  dimensionless  form: 


with  the  boundary  conditions:  (see  Fig.  1) 

p  (e  =  5,  n)  -  i 

P  (0  =  5  +  Of,  T|)  =  1 

p  (e  =  e,  n  =  +  j)  -  i  [2] 

Equation  [l]  can  be  reduced,  by  finite  difference  techniques,  to  the 
difference  equation: 


F?  .  -  2E,  .  P,  .  +  F.  .  —  0 


[3] 


Solving  equation  (3]  for  P.  .,  the  pressure  at  any  point,  yields 

J 

.  "  Ei,j  + 

where  the  coefficients,  E,  ,  and  F.  .  can  be  expressed  as: 

Ei,  j  =  2D(A0)2  EPi+l,J  +  Pi-1,  J  +  4D(A9)Hi  (S)±  [Pi+1, j  "  Pi-l,j] 
+  S)(to2^  [Pi, J+l  *  Pi, j-1  ] 


U] 


[5] 
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—  - [p  -  p  ]  — 1 —  rp  _  p  -|2 

2D(A6)H2  i+1>J  We)2  Lli+l,j  Pi-l,JJ 

-  - ^(gf  [P±  1+,  -  P±  ,  J2 

4D(At))2'L'  ijJ+l 


where 


[6] 


D  =  2 


[7] 


In  an  effort  to  facilitate  the  programming,  these  equations  were  rearranged 
into  the  following  form: 


Fi,J  "  CC  ^RiAi,J 

-  A2  -  TB2  ] 

[8] 

ei,j  -  °°  tcL + 1  si,/ 

‘i,j  +  -  5  Bisi] 

[9] 

where 

A 

Pi+1..1  "  pi-l..i 

R  A 

2 

Bi,j 

P  -  P 

i.,1+1 _ l.-i-l 

2 

s 

i  Ht\  d94 

* 

Ci,J 

P  +  P 

-ltlx.1 _ -1-1. .1 

2 

*i  -(i)2(^)2 

n 

P  +  P 

1-1+1  Pi.1-1 

C10] 

2 

CC  2(1  +  MN) 
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Equation  [4]  applied  to  all  points  of  the  bearing  grid  yields 
a  set  of  simultaneous  equations  which  can  be  solved  by  an  iterative  pro¬ 
cess.  In  the  preceding  equations,  the  pressure  at  any  point  is  defined 
as  a  function  of  the  pressures  at  the  four  adjacent  points.  An  arbitrary 
starting  value  is  chosen  for  the  pressure  at  all  grid  points  and  the 
process  is  repeated  until  the  sequential  adjustment  of  the  pressure  at 
each  iteration  is  smaller  than  a  pre-established  amount. 

At  this  point  it  is  important  to  make  a  basic  decision  between 
two  methods  of  attack: 

a)  Fix  the  geometrical  characteristics  of  the  bearing 
film  and  the  relative  speed  and  let  the  pressure 
distribution  and  center  of  pressure  be  a  consequence. 

b)  Fix  the  geometrical  characteristics  of  the  shaft 
and  bearing  and  the  position  of  the  center  of 
pressure  and  let  the  bearing  be  free  to  assume  any 
tilt  as  to  comply  with  the  specified  data. 

Method  a)  was  chosen  because  it  is  much  simpler  to  program, 
converges  more  easily,  does  not  get  involved  with  double-valued  problems, 
and  usually  gives  an  answer  in  half  the  time  employed  by  method  b).  The 
adopted  method  has  the  disadvantage  of  not  producing  direct  data  for  a 
specified  center  of  pressure  position  as  would  be  desirable  for  the 
solution  of  a  single  pivot ed-pad  bearing  problem.  However,  since  the 
production  of  extensive  field  maps  was  the  aim  of  the  project,  the  data 
necessary  for  the  study  of  pivoted-pad  bearings  can  be  easily  obtained 
by  cross-plotting. 

If  the  production  of  a  limited  set  of  data  for  a  particular 
pivot  pad  geometry  is  desired,  it  is  more  convenient  to  program  the  pro¬ 
blem  according  to  method  b). 
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As  is  common  with  iterative  processes,  the  determination  and 
careful  control  of  the  rate  of  convergence  is  vital  to  successful  pro¬ 
duction  runs.  Moreover,  if  the  produced  data  must  be  employed  in 
extensive  cross-plotting,  smoothness  and  accuracy  are  essential.  The 
criterion  for  truncation  must  be  so  selected  as  to  stop  all  iteration 
processes  at  a  uniform  "distance"  from  the  exact  answer.  Since  the 
exact  answer  is  not  known,  fictitious  criteria  must  be  selected  which, 
nonetheless,  are  required  to  be  equally  effective. 

The  following  presentation  expounds  one  such  criterion  based 
on  the  exponential  nature  of  the  solutions  generated  by  this  iterative 
technique.  The  convergence  of  the  selected  iterative  process  is  determined 
in  the  following  manner: 

Let  the  "iteration  increment"  at  step  k  be  defined  as 


This  is,  in  essence,  the  average  growth  of  the  pressure  at  a  non-boundary 
point  within  the  bearing.  Let  the  "total  growth  indicator"  at  step  k 
be  defined  as 


bb, 


k 

=  £ 

x=i 


sa. 


[12] 


This  indicator  measures  the  average  total  growth  of  the  pressure  profile 
after  k  iterations. 

The  successive  iterations  of  the  solution  of  Reynolds'  equation 
as  produced  by  either  the  "simultaneous  displacement"  or  "overrelaxation" 
techniques  have  been  found  to  behave  locally  exponentially  with  respect 
to  the  step  number  k.  This  fact  can  be  used  for  the  individual  extrapolation 
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of  local  pressure  values  toward  the  asymptotic  solution  in  a  manner  akin 
to  the  "Atkin  process".  Moreover,  it  can  be  speculated  that,  if  the 
individual  pressure  points  behave  exponentially,  also  the  total  growth 
indicator  might  show  the  same  effect.  Then  the  following  formula  can  be 
fitted  to  the  data: 

bb£  +  Ae_B1  =  bb.  [13] 


where 

bb^  =  the  total  growth  indicator  at  the  Ath  iteration 

bb^  =*  the  total  growth  indicator  after  an  infinite 
number  of  stable  iterations 

A,B  =  constants 

A  -  iteration  number. 

Applying  Equation  [13]  to  three  sequential  steps  and  eliminating 
the  constants,  one  can  obtain  the  following  relation: 


[14] 


For  the  particular  case  where  the  steps  are  successive  or  equally  spaced 


then  the  right  hand  side  of  equation  [14]  becomes 


Applying  this  to  equation  [14]  we  have 

bbl  "  bb?  \?  bbi  “  bb, 

*^  +  1)  -HfraJ  +  i. 


[16] 


[17] 
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Solving  Equation  [17]  for  bb^  results  in 


(bj2  -  b1  b? 


2b, 


9 


[18] 


Note  that  the  subscripts  1,  2  and  3  in  equations  [14]  through  [18]  denote 
three  sequential  iterations  of  the  process,  not  necessarily  the  1st, 

2nd  and  3rd  iterations.  (Figure  2  shows  the  plots  of  the  extrapolated 
growth  indicator  and  the  total  growth  indicator  vs.  the  iteration  number 
for  a  typical  case). 


-rr- - - 1 


bb_  -  bb 


bb. 


*<6 


[19] 


where  6  is  a  number  produced  by  a  predetermined  convergence  criterion. 

The  advantage  afforded  by  this  technique  consist  of  the  fact 
that  bb^  -  bbj  is  a  measure  of  how  "far"  the  iteration  solution  is  from 
the  exact  one.  Of  course,  it  could  be  observed  that  the  solution  does 
not  exactly  behave  according  to  the  above-mentioned,  exponential  law  and 
that  bbOT  is  indeed  bbw  -  bb^  (k),  k  *  iteration  number.  However,  but 
for  the  first  few  steps,  bb,,,  is  a  very  weak  function  of  the  step  number 
and  converges  to  the  asymptote  from  above  as  opposed  to  bb  which  converges 
to  the  asymptote  from  below. 

The  great  advantage  of  this  truncation  method  is  that  all  pro¬ 
cesses  are  stopped  at  approximately  the  same  distance  from  the  exact 
answer.  In  contrast  with  this,  most  iteration  techniques  in  the  literature 
are  stopped  when  the  iteration  increment  is  below  a  given  limit.  This 
casues  iteration  truncation  errors  which  vary  from  case  to  case  according 
to  the  exponential  behavior  assumed  by  the  solution.  This  technique  was 
experimented  with  at  the  start  of  the  project  and  produced  results  which 
could  not  successfully  be  cross-plotted  except  for  the  most  regular  regions 
of  the  field  maps. 
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FIG.  2  TOTAL  GROWTH  INDICATOR  AND  EXTRAPOLATED  GROWTH  INDICATOR  vs. 

ITERATION  NUMBER 
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Numerical  Stability  and  Pressure  Extrapolation 

The  selection  of  an  iteration  method  was  made  along  the 
following  lines: 

(a)  Simultaneous  displacement  techniques  were  chosen 
over  successive  displacement  ones  due  to  poor 
properties  presented  by  the  latter  in  connection 
with  the  numerical  stability  program. 

(b)  A  process  which  could  be  suitable  for  both  over¬ 
relaxation  and  underrelaxation  had  to  be  selected 
due  to  the  uncertain  numerical  stability  behavior 
of  the  solution. 

Processes  such  as  "simultaneous  relaxation  implicit  by  lines" 

and  the  "alternating  gradient"  method  (same  as  the  relaxation  implicit 

by  lines  but  in  alternating  directions)  were  tried  on  a  full  scale  basis 

3 

but  produced  mixed  results.  On  that  basis  they  were  discarded. 

Then,  the  growth  of  the  pressure  at  a  particular  grid  point 
is  controlled  by  the  following  equation: 

-  vfij1'  +  -  *>  ko:1 

_(m+l) 

where  P^  j  is  the  assigned  value  of  the  pressure  at  the  point  i,j  for 
the  (m+l5th  iteration.  p(®t^  is  the  calculated  value  of  the  pressure 

Ip  J  -j. 

at  the  point  i,  j  for  the  (m+l)th  iteration,  P7  .  is  the  value  of  the 

/  \  J 

pressure  at  the  point  i,J  for  the  (m)th  iteration. 

From  equation  [20 J,  it  can  be  seen  that  for  y  ■  1,  the  assigned 
value  of  the  pressure  at  the  point  i,J  is  identical  to  the  calculated 
value. 


Numerical  instability  in  the  iterative  process  is  detected  by 
an  increase  of  the  value  of  the  iteration  increment  aa^  (as  defined  in 
equation  [ll]).  This  term,  normally  behaves  as  a  exponential  and  only 
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when  numerical  instability  is  present  does  it  increase  rather  than  decrease. 
Usually,  instability  can  be  eliminated  by  lowering  sufficiently  the  value 
of  y.  This  attenuates  the  effects  of  erratic  higher  harmonics  of  the 
solution  and  induces  stability.  As  a  result  of  lowering  the  value  of 
y,  the  iteration  proceeds  at  a  slower  but  more  stable  rate.  Obviously 
the  speed  of  convergence  is  greatest  when  the  highest  value  of  y  for 
stable  iteration  is  used.  This  statement  holds  true  in  relaxation  processes 
for  most  elliptic  equations  departing  from  the  simple  la  place  equation. 
Indeed,  for  ^  P  -  0,  throughout  the  extensive  region  of  numerical  stability 
there  exists  a  value  of  y  which  provides  the  greatest  rate  of  convergence. 
This  value  of  y  is  called  "optimum  overrelaxation  factor"  and  the  process 
"optimum  overrelaxation".  However,  in  equations  noticeably  departing 
from  the  simple  Laplace  equation  the  stability  problem  dominates  and 
optimum  overrelaxation  cannot  be  achieved.  Indeed  it  is  only  in  rare 
cases  that  any  value  of  y  greater  than  one  can  be  used. 

With  all  of  the  previously  mentioned  conditions  imposed,  the 
number  of  iterations  required  for  the  solution  of  the  system  of  Reynolds' 
difference  equations  varies  depending  upon  the  bearing  parameters,  the 
number  of  grid  points,  and  the  numerical  stability  of  the  iterative  pro¬ 
cess.  In  order  to  reduce  the  number  of  required  iterations,  the  exponential 
behavior  of  the  solution  was  exploited.  The  use  of  the  exponential 
extrapolation  is  successful  in  most  cases  but  always  had  to  be  applied 
on  a  trial  basis  due  to  the  occasional  accentuation  of  local  irregularities 
into  large  errors.  Local  exponential  extrapolation  was  accomplished  by 
the  following  equation: 


p(n) 


[21] 


where 


p(nj!  is  the  pressure  of  the  point,  i,j  at  the  nth 
iteration, 


is  the  exponential  asymptote  of  the  pressure 
of  the  point  i,j 


-  14  - 


THE  FRANKLIN  INSTITUTE  •  Laboratories  for  Research  and  Development 

I-A2049-18 


A,B  are  constants. 


t  is  the  time  at  the  nth  iteration, 
n 

Use  of  the  same  technique  that  was  applied  to  equations  [14]  through  [18] 
results  in  the  following  expression: 


P(“) 


rp(2 ) 'j  p(l)  p(3) 


C22] 


The  superscripts  (l),  (2)  and  (3)  denote  three  equally  spaced  iterations 

in  this  sequence  at  any  point  in  the  iterative  process  and  not  necessarily 

00 

the  first,  second  and  third  ones.  Usually  P.  .  is  only  an  approximation 
to  the  exact  solution  of  the  governing  difference  equation  and  more 
relaxation  steps  have  to  be  applied  in  order  to  satisfy  the  desired 
truncation  criterion. 

\  . 

Judicious  application  of  the  extrapolation  routine  has  been 
found  to  be  the  most  useful  single  expedient  in  the  solution  of  this  class 
of  problems.  Extreme  care  must  be  exercised  in  extrapolating  only  grids 
which  are  proceeding  very  smoothly  and  which  do  not  present  any  irregularities 
or  "seeds  of  instability".  In  view  of  what  has  been  said  above,  it  can 
be  concluded  that  the  most  convenient  technique  to  follow  is  to  adjust 
grid  size  and  relaxation  factor  (y)  so  as  to  start  the  solution  on  a 
very  stable  basis  (values  of  y  corresponding  to  very  underrelaxed  conditions). 
A  few  iterations  will  then  set  the  basis  for  extremely  successful  extra¬ 
polations.  The  stable  relaxation  process  can  be  applied  again,  followed 
after  a  few  steps  by  another  extrapolation,  and  so  forth.  The  establish¬ 
ment  of  a  proper  sequence  is  rather  empirical  but  important  for  economization 
of  valuable  machine  time  in  large  production  runs. 


It  is  important  to  note  that  experience  with  extrapolation  was 
less  successful  with  geometries  which  involved  convergent-divergent  films 
than  with  the  simply  convergent  or  simply  divergent  ones.  An  example  of 
the  use  of  pressure  extrapolation  can  be  3een  in  Figure  2.  At  the  25th 
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iteration,  the  extrapolation  routine  was  applied  and  the  step  in  the 
total  growth  indicator  curve  is  indicative  of  the  reduction  of  the  number 
of  required  iterations. 


Steady-State  Bearing  Characteristics 


Once  satisfactory  convergence  of  the  pressure  profile  over  the 
bearing  surface  is  obtained,  it  is  possible  to  solve  for  the  bearing 
characteristics.  The  load  capacity  for  the  bearing  is  defined  as 


where 


[233 


[24] 


[25] 


The  frictional  coefficient  is  defined  as 


These  integrations  are  performed  by  Simpson's  rule. 

Moreover,  the  knowledge  of  the  load  components  can  be  utilized 
to  establish  the  position  of  the  center  of  pressure  with  respect  to  the 
bearing. 
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3.  GENERAL  CONSIDERATIONS  FOR  APPLICATION  OF  COMPUTER  DATA 

The  characteristics  of  gas-lubricated,  self-acting,  cylindrical, 
partial-arc  journal  bearings  have  wherever  possible  bean  identified  by 
standard  terminology.  However,  there  are  certain  parameters  used  in  this 
report  that  are  unique  to  these  bearings.  These  parameters  are  defined 
in  the  text  that  follows: 

The  Pivot  Circle  Clearance.  C' 

This  parameter  establishes  the  relative  position  of  the  pivotal 
point  (taken  on  the  bearing  surface)  of  tilting  pad  shoes  with  respect 
to  the  shaft.  We  define 


“  Rpivot  ‘  Rshaft  k’ 

where  R^^  is  the  distance  from  the  pivotal  point  to  a  properly  selected 
point  fixed  with  respect  to  the  bearing  mounts.  For  the  case  of  three 
shoe  arrangements  R^^  is  the  radius  of  the  circle  passing  through  the 
three  pivot  points,  ^8haft  is  radius. 

The  Ground  Clearance  C 

Is  defined  as  the  difference  in  radii  of  the  bearing  and  shaft 

surfaces 

C  ”  bearing  "  Rshaft 

The  Dimensionless  Pivotal  Film  Thickness  H 

—  -  -  ■  - -  p 

The  film  thickness  at  the  center  of  pressure  or  pivot  position 
is  an  important  parameter  in  determining  the  operation  of  the  bearing. 

For  a  fixed  shaft  position  it  will  remain  a  constant  regardless  of  the 
attitude  of  the  bearing  pad  with  respect  to  the  pivotal  circle.  In 
agreement  with  the  geometry  of  Figure  1  this  film  thickness  can  be 
written  as: 
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Hp  ■  (C  +  e  cos  +  *})/C  [28] 

The  Pad  Position  Angle  S 

5  is  defined  as  the  angle  drawn  from  the  line  of  centers  of 
the  bearing  and  shaft  circles  to  the  leading  edge  of  the  pad. 

Trailing  Edge  Bearing  Number,  A^, 

The  parameter  "trailing  edge  bearing  number"  is  defined  alge¬ 
braically  as 

where 

hj,  **  C[l  +  c  cos  (5  +  »)]  [30] 


[31] 


This  number  more  closely  represents  the  operation  of  an  individual  shoe 
from  the  conventional  compressibility  bearing  number  point  of  view. 

The  Pivotal  Bearing  Number  Ap 

This  form  of  the  bearing  number  is  based  on  the  clearance  that 
exists  at  the  center  of  pressure,  or  pivot.  Algebraically,  it  can  be 
expressed  as 

A  =  (±f  [32] 

P  Pa  \h?)  ’ 


where 

hp  —  C[1  +  c  cos  (§  +  0)] 


[33  3 


[34] 
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The  Clearance  Ratio  C'/C 

This  parameter  is  defined  as  the  ratio  of  the  pivotal  circle 
clearance  to  the  ground  clearance.  It  is  important  because  it  measures 
the  amount  of  preload  existing  due  to  the  initial  pivot  circle  adjust¬ 
ment.  As  an  example,  consider  a  three  shoe  bearing  system  in  which  each 
shoe  has  an  angular  extent  of  120°.  The  film  thickness  at  the  pivots 
can  be  expressed  as 

hPl  *  C'  Cl  +  «'  cos  X]  -  hpi  =  CHpi  *  [35] 

hp2  -  C«  [1  +  «'  cos  (X  +  120°)]=  hp2  =  CHpg  [36] 

hp3  =  C'  [1  +  «'  cos  (X  +  240°)]  =*  hpj  =  CHp3  [37] 

Sunming  Equations  [35 3,  [36]  and  [373  yields  the  following: 

npi  +  hP2  +  ^3  "  c7  *3  +  s'fcos  X  +  coe(x  +  120 +  cos(x  +  240  [38] 


Performing  the  trigonometric  substitutions  yields 

Cl  _  HP1  +  HP2  +  ^3 
C  3 


[39] 


This  formula  cannot  be  used  if  the  spacing  between  pivots  is  not  equal. 
Indeed,  for  the  more  general  case  in  which  symmetry  exists  about  a 
vertical  plane  only,  formula  [293  becomes 


c,  Hpi  +  Hp2  +  (2  cos  P)Hp3 
C  **  *  2(1  +  cos  0) 


[39a] 


where  2@  is  the  angle  between  pivots  1  and  2. 

Pivotal  Circle  Eccentricity  Ratio  «' 

The  pivotal  circle  eccentricity  ratio  is  defined  as  the  relative 
location  of  the  shaft  center  with  reference  to  the  pivotal  circle  center. 

•a 

X  -  5  +  0 
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4.  DESIGN  DATA 

The  data  obtained  from  this  analysis  can  be  presented  in  many 
different  arrangements.  Three  basic  types  of  plots  are  presented  in  this 
report.  They  consist  of  plots  of  (a)  load  coefficient  (C^)  vs.  center 
of  pressure  position  (d/o)j  (h)  pivotal  film  thickness  (Hp)  vs.  center 
of  pressure  position  (#/«);  (c)  trailing  film  thickness  (H^)  vs.  center 
of  pressure  position  (®/a). 

All  plots  of  the  (a),  (b)  and  (c)  variety  contain  lines  of 
constant  c  and  lines  of  constant  5.  Since 

CL  *  f(t,  5,  A,  R/L,  or) 

and 

4>/a  -  g(e,  5,  A,  R/L,  of) 

etc.,  the  best  a  two  dimensional  plot  can  do  is  to  present  vs.  on 
lines  of  constant  s  and  lines  of  constant  ?  while  A,  R/L,  a  must  have 
assigned  values  on  each  graph.  It  should  be  noticed  that  for  compressible 
lubrication  the  speed  parameter  A  and  the  load  parameter  cannot  be 
combined  in  a  manner  similar  to  the  Sommerfeld  number  of  incompressible 
lubrication.  This  fact  is  due  to  the  non-linearity  of  the  governing 
equation. 

It  is  clear  that  the  amount  of  information  necessary  to  exhaust 
all  possible  interesting  geometries  and  running  conditions  can  assume 
colossal  proportions.  Our  attitude  was  to  investigate  rather  thoroughly 
the  geometries  which  were  closely  connected  with  other  vx>rk  at  The 
Franklin  Institute  and  also  to  make  the  program  available  to  anyone 
interested  in  obtaining  more  data.  Two  values  of  angular  extent  were 
investigated:  a  =  120°  and  a  =  94.5°.  The  first  angle  closely  approximates 
symmetrical  three-shoe  arrangements.  The  second  was  actually  used  for  a 
practical  design  where  the  unidirectionality  of  the  load  made  it  preferable 
to  space  the  lower  pivots  100  degrees  apart  rather  than  at  an  angle  of  120°. 
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For  the  a  =  120°  runs  the  aspect  ratio  parameter  R/L  was  given 
values  of  1/3  and  1/2*  For  the  a  m  94-5°  case  the  values  of  R/L  employed 
were  the  ones  dictated  by  the  actual  parts  of  the  experimental  machine 
under  test  at  The  Franklin  Institute  (AEC  Contract  No.  AT( 30-1 )-2 512, 

Task  3)j  these  were  R/L  -  0.6061  and  R/L  *»  0.4041  (corresponding  to  width 
to  length  ratios  of  the  developed  surface  of  1.0  and  1.5  respectively). 

This  report  contains  the  part  of  the  results  which  were  most 
completely  organized  up  to  the  present  time:  the  94.5°  shoes  with 
R/L  =  .6061.  The  actual  tabulation  of  data  is  very  lengthy  and  is  not 
contained  herein.  However,  it  is  worth  mentioning  that  all  computer 
data  plotted  so  smoothly  that  no  point  would  fall  off  the  lines  presented 
in  the  figures.  Another  important  factor  in  assessing  the  value  of  the 
present  solution  is  that  actual  checks  of  theory  against  experiment  were 
performed  with  extremely  favorable  results.  This  part  of  the  work 
together  with  a  thorough  analysis  of  the  data  and  evaluations  of  design 
techniques  will  be  presented  in  a  subsequent  report. 

Figures  3  through  17  contain  the  above  mentioned  field  maps  for 
values  of  the  bearing  parameter  A  of  1.5,  2.5,  3.0,  3.5  and  4.0.  These 
values  of  A  are  representative  of  typical  operating  conditions.  The 
effects  of  compressibility  present  in  a  bearing  geometry  are  not  com¬ 
pletely  represented  by  the  value  of  A.  This  is  due  to  the  fact  that  A 
is  based  on  the  ground  clearance  "C"  between  the  bearing  pad  and  the 
shaft.  According  to  the  values  of  e  and  5  the  value  of  the  actual 
clearance  and  particularly  the  value  of  the  trailing  edge  clearance  can 
be  much  lower.  For  this  reason  the  value  of  A^  (the  compressibility 
parameter  based  on  the  trailing  edge  film  thickness)  more  closely  represent 
the  extent  of  the  compressibility  effects  in  the  bearing  film. 
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FIG.  3  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL  BEARING 
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FIG.  4  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL  BEARING 

FOR  A  =  1.5 
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FIG.  7  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL  BEARING 

FOR  A  =  2.5 
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FIG.  9  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL  BEARING 

FOR  A  =*  3.0 


£ 
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FIG.  10  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL  BEARING 

FOR  A  =  3 .0 
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FIG.  12  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL  HEARING 

FOR  A  =  3.5 
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FIG.  15  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL  BEARING 

FOR  A  *  4.0 
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FIG.  16  THEORETICAL  PERFORMANCE  OF  PIVOTED  PARTIAL  GAS  JOURNAL 
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5.  USE  OF  FIELD  MAPS 

The  characteristics  of  various  bearing  arrangements  can  be 
predicted  by  proper  use  of  Charts  3  through  17.  In  all  the  calculations 
frictional  forces  will  be  neglected  in  comparison  to  pressure  forces. 

All  fixed  geometries  can  be  easily  analyzed  because,  for  any  position 
of  the  shaft  within  the  bearing,  the  exact  geometry  of  each  bearing  film 
is  known  (*  and  5).  Then  charts  3,  6,  9,  12  and  15  give  directly  the 
loads  and  points  of  application  of  the  loads  in  each  bearing  segment. 
Vector  addition  of  the  loads  gives  the  total  load  and  attitude  angle. 

In  particular,  axial  groove  cylindrical  bearings  have  pads 
characterized  by  a  conmon  value  of  e  equal  to  the  eccentricity  ratio  of 
the  shaft  within  the  bearing  circle.  On  a  constant  s  line  the  pads 
correspond  to  values  of  5  starting  from  an  arbitrary  number  and 
spaced  at  intervals  equal  to  a  (angular  extent  of  the  pads).  Vector 
addition  of  the  loads  will  produce  a  total  load  and  attitude  angle 
corresponding  to  each  selection  of 

CL  -  f(So> 

6  =  attitude  angle  =  g(5Q) 

In  general,  it  will  be  of  interest  to  select  only  the  solution  for  which 

T 

SQ  -  6  =  50  =  angle  between  load  line  and  leading  edge  of  first  pad 
(see  Fig.  18).  The  independence  of  A  from  the  shaft  position  is  a  great 
asset  in  these  calculations.  Indeed  bearing  parameters  based  on  any 
of  the  individual  pad  clearances  would  require  cross-plotting  of  the  data 
to  have  A  on  one  of  the  coordinate  axes  and  maps  of  constant  eccentricity. 
Supplementary  nomographs  relating  A  to  geometry  would  also  be  needed. 
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In  the  treatment  of  pivot ed-pad  bearinge  the  geometry  becomee 
more  complicated.  Figure  19  shows  this  type  of  bearing  arrangement  for 
the  three  shoe  case  with  symmetry  about  a  vertical  diameter. 

The  invar lent  characteristic  of  a  working  pivoted  pad  ie  the 
position  of  its  pivot  and-at  any  journal  speed  -  the  value  of  A.  This 
means  that  all  information  possibly  pertaining  to  the  problem  is  contained 
on  appropriate  vertical  lines  of  Figures  3  through  17.  Then,  for  any 
position  of  the  shaft  within  the  pivot  circle  (characterized  by  the 
eccentricity  ratio  s’  and  the  attitude  angle  6’),  the  pivotal  clearances 
Hp  are  known.  One  of  the  graphs  of  Hp  vs.  0/a  (the  one  for  the  appropriate 
value  of  A)  will  translate  this  information  into  the  values  of  e  and  5 
for  each  pad.  The  graphs  of  vs.  <t>/a  will  now  be  used  to  obtain  the 
vector  loads  through  each  of  the  pivots.  Thd  vectorial  addition  of  the 
component  loads  produces  the  total  load  vector.  For  any  eccentricity 
ratio  e',  trial-and-error  techniques  have  to  be  used  in  selecting  a  value 
of  6'  resulting  in  a  total  load  in  the  wanted  direction. 

Considerable  simplification  is  afforded  by  special  geometries 
such  as  the  one  with  symmetry  about  a  vertical  diameter.  In  this  case 
the  attitude  angle  is  known  to  be  zero  and  the  trial-and-error  procedure 
is  no  longer  necessary.  As  an  example  let  us  consider  the  problem  of 
evaluating  the  total  vs.  *'  relation  at  any  particular  value  of  the 
bearing  3peed  parameter  A.  For  this  case  let  us  use,  (see  Fig.  19) 

a  =  94.5°  (Three  equal  pads) 

R/L  =  0.606l 
A  =  3.5 
P  -  50° 

Then  the  pivot  clearances  are  given  by 

Kpi  -  np2  -  ^  [!  +  «'  COS  (tt  -  P)]  *  2jr  Cl  +  0.643  *’] 
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JOURNAL  BEARING 


FIG.  19  PIVOTED  PAD  JOURNAL  CONFIGURATION  WITH  THREE  SHOES 
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^3  “  c“  t1+*'C08  °°3  -  y  Ci  +  f  ] 


From  the  values  of  C^,  C^,  C^y  the  total  load  carrying  capacity  is 
obtained  as 

°LT  "  20U  “•  "  ‘  CL3 


CLT  -  1-2S6  CU  -  °L3 

The  values  of  G^,  C^,  G^  are  obtained  by  (a)  entering 
Figure  13  with  the  proper  values  of  Hp  and  0/a,  reading  out  the  corres¬ 
ponding  * ' s  and  ?'a  and  (b)  entering  Figure  12  with  0/a,  c  and  £  and 
read  out  C^.  After  step  (a)  the  values  of  t,  5,  0/a  can  be  used  to  enter 
Graph  14  and  read  out  H_.  These  calculations  are  shown  in  Table  I  for 

c f  1 

values  of  -g-  equal  to  0.6,  0.8  and  1.0  and  the  results  are  plotted  in 
Figure  20. 

An  example  of  the  use  of  a  plot  such  as  Figure  20  is  given  by 
the  following  problem: 

Given:  Rotor  weight  =  80  lb 

Rotor  radius  =  2  in. 

Pad  axial  length  =  3*3  in. 

C  *  1.5  x  10"3  in. 

Angular  velocity  =  17,650  rpm 

Find:  Trailing  edge  film  thickness,  running  eccentricity, 

and  stiffness  corresponding  to  various  values  of  pivot  circle  setting. 


Procedure : 


c  =  _ _ 42 - -  o 

L  P&RL  14.7  x  2  x  3.3 
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A  -  &asi  ,  $  *  (2.61  S  IQ'9)  XJL7.W6Q  X  2n(2)2 
Pac2  14.7  x  (1.5  x  10"3)2 


A  -  3.5 

Using  Figure  20 


C' 

c 

«'  . 

Sl. 

_h_ 

(in) 

acT 

_ 4 

a«' 

deW 

Per  Bearing 

0.6 

0.26 

0.379 

0.569 

x  ioL3.. 

1.14 

1.23  x  105 

0.8 

0.40 

0.427 

0.640 

x  10"3 

0.91 

0.736  x  10? 

1.0 

0.54 

0.460 

0.690 

x  10'3 

0.80 

0.517  x  105 

where  the  stiffness  (per  bearing)  is 


aw 

do 


a  yi 

at*  (_c  x  (C'/cLl  ai'L 


From  the  results  it  can  be  seen  that  the  preload  parameter 
C'/C  has  a  marked  effect  on  the  operation  of  pivot ed-pad  gas  bearings. 
Particularly  noticeable  is  the  effect  of  C'/C  on  the  stiffness  which  can 
be  interpreted  as  a  valuable  tool  for  aposteriori  adjustment  of  the 
rotor  critical  speed.  Since  it  is  well  known  that  the  self-eoccited 
whirl  threshold  speed  is  some  multiple  of  the  rotor  critical  speed 
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(approximately  twice  for  email  attitude  angles ) ,  the  whirl  instability 
region  can  be  modified  within  limits  by  proper  setting  of  the  pivot 
circle  clearance. 


Dudley  D.  Puller 
Project  Engineer 


Approved  by: 
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